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Hi and CO observations of Arp 104: a spiral-elliptical interacting 
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ABSTRACT 

We present data probing the spatial and kinematical distribution of both the atomic (H I) and 
molecular (CO) gas in NGC 5218, the late-type barred spiral galaxy in the spiral-elliptical 
interacting pair, Arp 104. We consider these data in conjunction with far-infrared and radio 
continuum data, and A^-body simulations, to study the galaxies interactions, and the star for- 
mation properties of NGC 5218. We use these data to assess the importance of the bar and 
tidal interaction on the evolution of NGC 5218, and the extent to which the tidal interaction 
may have been important in triggering the bar. The molecular gas distribution of NGC 5218 
appears to have been strongly affected by the bar; the distribution is centrally condensed with 
a very large surface density in the central region. The A^-body simulations indicate a timescale 
since perigalacticon of ~ 3 x 10 8 yr, which is consistent with the interaction having triggered 
or enhanced the bar potential in NGC 5218, leading to inflow and the large central molecular 
gas density observed. Whilst NGC 5218 appears to be undergoing active star formation, its 
star formation efficiency is comparable to a 'normal' SBb galaxy. We propose that this system 
may be on the brink of a more active phase of star formation. 
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1 INTRODUCTION 

Barred galaxies represent an important sub-sample of systems in 
any study of galaxy interactions and induced star formation. Bar 
formation is attributed to one of two possible mechanisms; bars 
can form spontaneously as a result of an intrinsic gravitational 
instability, or they can be produced by tidal disturbance of the 
galactic potential during an interaction. Both mechanisms are sup- 
ported by numerical simulations. Massive stellar disks with small 
velocity dispersions are intrinsically unstable to bar formation 
and will develop a strong bar structure without external perturba- 
tion w ithin a few rotation periods (e.g. jAthanassoula & Sellwooj 
1986). Similarly, external perturbations can induc e bar structures 
in galactic disks that are stable in iso l ation (e.g., Noguchil 
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Geri n. Combes & Athanassoulal Il990l : iMihos & Hernquisl 
hereafter MH96) 



1987; 



1996, 



Observationally, iKeel et all dl996h used K-band images to 
study the bar fraction in a sample of pairs, fin ding 80% of the pairs 
to disp lay signs of bar formation. lElmegreen. Elmegreen & Bellini 
( 1990) also find a high fraction of barred galaxies in binary systems 
(50 ±5%), higher than in the field (32 ±9% ) or group environments 
(28 ±2%). However. Ivan den Berghl fe002h . using a sample of 930 
galaxies found no evidence for a dependence of bar frequency on 
galaxy environment, consistent with the nature of the host galaxy 



being the most important factor in bar formation. Taken together 
these results suggest that bar formation is strongly dependent on 
the host galaxy, but that pair-wise interactions do induce bars in 
many systems. 

Galaxy bars represent a principal mode of disruption to the 
axisymmetric potential of spiral galaxies, providing a mec hanism 
of gas inflow toward the central regions. The simulations of Mihos 
( 1994) and MH96 indicate that bars can play a particularly impor- 
tant role in fly-by encounters, where they provide a dominant trig- 
ger for central gas concentration and ensuing star formation. In ad- 
dition, Jhejvigjrrjroportonof^alaxi*^ 

bars dEskridge et al. 2000; Mene ndez-Delmestre et al.l2 004) means 
that the response of a barred potential to an external perturbation is 
also of interest. The relationship between bars and interactions and 
the evolution of the ISM in barred systems is crucial in attaining a 
better understanding of the star formation response of interacting 
galaxies. 

B oth theoretical (MH 96) and observational ^Sakamoto et al.l 

Il999al ; ISheth et al.l [2005) evidence exists pointing to enhanced 
molecular gas concentration in barred systems (at least in early- 
type spirals). However, the relationship between bars and star for- 
mation remains less clear. Whilst bars do not appear to affect the 
disk-wide star-formation properties of galaxies, they can be impor- 
tant in triggering enhanced star formation in the central regions (for 
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a review see lKennicutd[l994l) . Star formation in barred galaxies 
manifests itself in a div erse range of morphologies (Philli pj 19961 ; 
ISheth et alj200cll2002h . A number of studies point to an enhanced 
level of star formation in the central region of early-type barred 
systems, altho ugh similar enhancements are not observed in late- 
type systems ( Huang et al.|[l996l iHo et alj[l997h . This distinction 
is consistent with molecular gas observations. Bars may well play 
a particularly importa nt role in the e volution of the most actively 
star-forming systems. IPhillipsI Il996h examined the Lpjr/Lb ratio 
in a sample of IR-bright galaxies and only detected a difference be- 
tween the barred and non-b arred systems for the mo st actively star- 
forming galaxies. Similarly, Hunt & Malkan ( 1999) f ound a highe r 
incidence of bars in starburst galaxies, and|fjawarden et al.l dl986h . 
who looked at the properties of intrinsically luminous galaxies, de- 
tected a mean IR luminosity three times larger in barred systems 
than their unbarred counterparts. 

In this paper we examine the spatial and kinematic dis- 
tribution of both the atomic and molecular gas associated with 
NGC 5218, the late-type barred spiral galaxy in the spiral- 
elliptical interacting pair Arp 104. This is one of a small 
sample of nine spiral-elliptical interacting pairs, which have 
been studied by Cullen, Al exander & Clemens! J20031 120061) and 
ICullen & Alex ander d2005l) . Each pair has an H I flux of at least 
5 Jy km s~ 1 , and Arp 104 has a relatively large separation between 
the galaxies, implying they are in the early state of interaction. We 
use these data in conjunction with yV-body simulations and data 
probing the star formation properties of this system to assess the 
impact of the bar on NGC 5218's evolution, and look at the extent 
to which the tidal interaction may have been important in trigger- 
ing its bar. In Section|2]we present details of our H I and CO obser- 
vations. In Section [3] we discuss the star-formation properties and 
efficiency of NGC 5218 and present the results of the simulations. 
In Section|4]we examine the importance of the bar in this galaxy. In 
Section [5] we discuss our results, looking first at the importance of 
the bar in terms of NGC 5218's evolution and subsequently at the 
likelihood that the bar has been triggered by the tidal interaction. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Optical morphology of Arp 104 

Arp 104, shown in Fig.[T] consists of two interacting galaxies, see 
NEE0: the southern galaxy, NGC 5216, is an E0 pec (2'.5 x lf5) and 
the northern galaxy, NGC 5218, is an SB(s)b? pec (1 '.8 x l'.3). At an 
assumed distance of 4 1 . 1 Mpc (from a recessional velocity of 2880 
km s , and using a Hubble constant of 70 km s Mpc -1 ), the 
projected major axis diameter (D25) ls 29.9 kpc (2.5) for NGC 5216 
and 21.5 kpc (L8) for NG C 5218. The projected separation of 
the tw o galaxies is 49.1 kpc. lSchombert. Wallin & Struck-Marcelll 
4l99Ch have conducted a multicolour photometric study of the op- 
tical tidal features of a sample of interacting galaxies, including 
Arp 104. They found the tidal features of Arp 104 are bluer than 
the central regions of both NGC 5218 and NGC 5216, but similar 
to the colours of the outer disk in NGC 5218. 



2.2 VLA 21-cm observations 

The H I observations of Arp 104 were made with the VLA in C- 
array configuration on 1994 December 24 and in D-array on the 
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1995 March 22. Both datasets were retrieved from the VLA archive 
and analysed. For both the C- and D-array observations, a correla- 
tor mode with online Hanning smoothing giving 63 channels with 
a velocity width per channel of 10.3 km s (48.8 kHz) was used. 
Both sets of observations were centred on RA 13 h 32 m 09?36, Dec 
+62° 44' 03f'3 (J2000), at a heliocentric radial velocity of 2900 
kms -1 . The total on-source integration times were 258 and 236 
minutes for the C and D-array observations respectively. In both 
cases on-source scans were interleaved with those of the phase cal- 
ibrator, 1400+621. For the C-array observations both 3C48 and 
3C286 were observed as flux calibrators during the observations, 
with assumed flux densities of 16.0 and 14.8 Jy respectively. For 
the D-array observations 3C286 was used as the flux calibrator with 
an assumed flux density of 14.8 Jy. 

Data reduction followed standard procedures using the Na- 
tional Radio Astronomy Observatory (NRAO) Astronomical Image 
Processing System (AIPS). After flagging, the data were binned 
to velocity resolution of 20.7 kms -1 to improve the signal to 
noise. Subtraction of the average of the line-free channels from 
all channels was used to remove the continuum emission. AIPS 
task IMAGR was used to produce naturally- weighted channel maps 
with an angular resolution 21"3 x 16l'4 for the C-array data and 
63 "9 x 54"4 for the D-array data. 



2.3 GMRT H I absorption observations 

In addition to the H I images from the VLA observations, 
NGC 5218 was also observed with the Giant Metrewave Radio 
Telescope (GMRT) on 2003 October 16 to obtain an Hi absorp- 
tion spectra towards its nuclear continuum emission. The GMRT 
is an apertu re synthesis telescope comprising of thirty 45-m dishes 
jRaol2002). The full width half maximum of the primary beam of 
the telescope at 1.4 GHz is ~ 25'. The backend setup employed 
was a 30 station FX correlator with 128 spectral channels over a 
4 MHz bandwidth (857 km s _1 ), giving a resolution of 31.25 kHz. 
The observations were centred at a heliocentric radial velocity of 
2880 km s . The total integration time on source was 174 min- 
utes. On-source scans were interleaved with those of the phase cal- 
ibrator, 1313+675. 3C286 was used as a flux calibrator with an 
assumed flux density of 14.8 Jy. The pointing centre was the core 
of NGC 5218: RA 13 h 32 m 10?50, Dec +62° 46' 04"0 (J2000). 
These data were binned to a velocity resolution of 40 km s and 
channel maps made using natural weighting, giving a spatial res- 
olution of 7"72 x 4"60. The flux density of the continuum source 
is approximately 2 5 mjy, consistent with the FIRST observation 
(Whi te etaljl997h . 



2.4 12 CO (J = 1 0) observations 

Observations of NGC 5218 in the 12 CO (J = 1-0) line at 
115.27 GHz were made using the 10-element Berkeley-Illinois- 
Maryland Association (BIMA) millimeter interferometer. Two us- 
able tracks were observed in C-configuration on 2003 April 28 and 
May 13. The BIMA correlator was configured to have a resolu- 
tion of 1.56 MHz (4 km s _1 ) over a total bandwidth of 368 MHz 
(959 kms -1 ). The primary beam at BIMA at this frequency is 
100" FWHM. 3C84 and 1638+573 were observed for passband 
and phase calibration respectively. Data calibration was performed 
in Miriad dSault. Teuben & Wrightll995h . Calibrated datasets were 
then read into AIPS, combined, averaged to a spectral resolution of 
20 km s and images synthesised. 
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Figure 1. Integrated H I column density from combined VLA C- and D-array data (Sfwhm = 27"4 x 24"6 at a position angle of 53?3) for Arp 104 overlayed 
on a digitised sky survey fi-band image. Contour levels are (0.05, 0.075, 0.1, 0.125, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4) Jy beam -1 km s" 1 . 



2.5 12 CO (7 = 2 1) observations 

Observations of the 12 CO (7 = 2-1) emission of NGC 5218 at 
230.54 GHz were made with the JCMT in 'service' mode on 2001 
November 30 and 2003 January 4. Observations were made on 
a 9 x 5 grid in RA and Dec, with 10" spacing, centred at RA 
13 h 32 m 11?8, Dec +62° 46' 05" (J2000). The 10" spacing of 
grid points produced a fully-sampled map over the area observed 
(JCMT beam 20"8 at 230.54 GHz). Sky subtraction was achieved 
by beam switching to a point 80" in declination from the reference 
position. Observations were centred on a heliocentric radial veloc- 
ity of 2880 km s -1 with a bandwidth of 920 MHz (1200 km s -1 ), 
split into channels separated by 625 kHz. The spectra were binned 
to a velocity resolution of 20.9 km s to match that of the VLA 
data. 

Standard calibration observations were made at the same time 
as the observations of NGC 5218. Data calibration followed the 
standard JCMT procedure, adopting the beam efficiency and for- 
ward scattering efficiency on the JCMT website yielding an T7f ss 
of 0.77. To convert between main beam temperature and flux it 
was assumed that the 12 CO (7 = 2—1) emission could be modelled 
as a point source. Whilst this is not strictly the case, it allows a 



lower limit to be placed on the 12 CO (7 = 2-1) flux for NGC 5218. 
The conversion between flux and main beam temperature, assum- 
ing a point source, is given by (5/Jy) = 1 8.4(7^ /K) (assuming 
the FWHM of the JCMT beam at 230 GHz is 20!' 8). 



3 RESULTS 

3.1 Star-formation rate 

The flux densities of NGC 5218 at 12, 25, 60 and 100 jlm, taken 
from the IRAS High Resolution Image Restoration Atlas l lCao et ail 
Il997t) . are 0.36, 0.92, 7.14 and 14.38 Jy res pectively. This yie lds 
a far-infrared (FIR) luminosity (40-122 um. lHelou etalll 19881) of 
8.37 x 10 36 W (2.2 x 10 10 L©). Using the calibration of lBelll ( l2003h 
we obtain a FIR star-formation rate (SFR) of 5.3 M© yr _1 , The 
infrared lumin osity (8-1000 tim) for NGC 5218 is 3.96 x 10 10 L© 
dPeraultlll987t) . 

Both NGC 5218 and NGC 5216 have d etected radio contin - 
uum emission recorded by the NVSS survey dCondon et "all 19981) . 
with 30.4 and 31.8 mjy respectively. The NVSS flux density of 
NGC 5218 corresponds to a lumi nosity of 6. 16x 10 21 W Hz -1 . 
Using the radio SFR calibration of iBelj ( 120031) this yields a SFR 
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Figure 2. Channel maps for the combined VLA C and D-array data (Qfwhm = 27"4 x 24" 6 at a position angle of 53?3) for Arp 104 overlayed on a digitised 
sky survey B-band image. Maps have been produced by averaging three adjacent channel maps together. Contour levels are every 0.005 Jy beam -1 km s~', 
omitting zero (negative contour are dashed and black, positive contours are solid and white). 
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Figure 3. The intensity weighted mean velocity field of H I emission from the combined VLA C- and D-array data (Sfwhm = 27"4 x 24"6 at a position angle 
of 53?3) for Arp 104. The colour scale goes from 2716 (blue) to 3093 km s" 1 (red) and contour levels are (2760, 2780, 2800 . . .3100) km s -1 . 



of 3.4 Mq yr , a factor of ~ 2 below the SFR obtained using the 
calibration of ICondon] d 1992b (7.6 Mq yr -1 ). For the purposes of 
this work we adopt the FIR estimated SFR of 5.3 M yr -1 . 



3.2 Atomic Hydrogen 

The combined VLA C- and D-array Hi map of Arp 104 (Fig.Q3 
reveals a tidal tail to the north-east of the late-type system 
(NGC 5218), and a tidal bridge extending southward toward the 
early-type companion (NGC 5216). Figure 1, like the other H I and 
CO results shown here, has not been corrected for the primary beam 
of the telescope, which is much larger than the scale of the emission 
that is seen. The tidal tail to the north-east forms a quarter-circle arc 
with a radius of approximately 4', corresponding to a physical ra- 
dius of 48 kpc. The tidal bridge to the south is approximately 5' 
long, corresponding to a physical length of 60 kpc. The bridge is 
oriented directly toward the elliptical, NGC 5216, with a small tail 
extending beyond the elliptical to the south-west. 

Using the lower resolution D-array data we detect an H I flux 
above a level of 0.035 Jy beam -1 km s -1 over an area of approx- 
imately 37 arcmin 2 . At a distance of 41 Mpc, this corresponds to 
a physical area of 5360 kpc 2 . The detected flux in this low resolu- 



tion map peaks at 0.98 Jy beam - km s , falling to around 0.05 
Jy beam -1 km s -1 in the more extended regions of the tidal tail. 

The total Hi flux of the Arp 104 system, obtained by adding 
all the flux in the low-resolution channel maps is 7.0 ± 0.5 Jy 
kms -1 , which is in good agreement with the single-dish flux, 
7.7 ±2.4 Jy kms -1 dTheureau et al.l ll998). The Hi mass, given 
by 

Mm — 2.36 x 10 5 £> 2 Sdv M , 

is (2.8 ±0.2) x 10 9 Mq. 

Approximately 30% of the H I emission detected in the D- 
array map appears to be associated with the late-type system, 
NGC 5218, with the remainder forming an extended tidal tail. 
This i s consistent with the trend found by Hibb ard & van Gorkoml 
dl996h . For a sample of four galaxies (Arp 295b, NGC 4676a, NGC 
4676b and NGC 520) considered to lie in the central region of the 
merger sequence, they find ~ 35% of the total atomic gas emission 
comes from the inner regions of the system; this percentage falls 
dramatically for more advanced mergers. 

Channel maps of the Hi emission are shown in Fig. [2] 
whilst Fig.|3]displays the intensity weighted mean velocity field of 
Arp 104. H I emission is detected in channels 7-25 corresponding 
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Figure 4. H I absorption spectra towards the central continuum source of 
NGC 5218. Data at a spatial resolution of 7"72 x 4"60 at a position angle 
of —40? 1 binned to a velocity resolution of 40 km s . 



to the velocity range 2716—3094 km s . The northern tidal tail has 
a velocity dispersion of ~ 270 km s with velocities in the range 
2821—3094 kms -1 . The tidal bridge to the south has a slightly 
smaller velocity dispersion of ~ 210 km s , with velocities in the 
range 2800—3010 km s . The largest velocity dispersion is ob- 
served for the gas still associated with NGC 5218. Here velocities 
range from 2716—3094 km s with the western side blueshifted. 

The central region of NGC 5218 is seen in absorption due 
to the presence of a bright radio continuum source in NGC 5218 
(see Fig.|4). The absorption profile is centred on a velocity of 2884 
km s . The detected H I emission believed to be still bound to the 
late-type system is distributed approximately symmetrically about 
this velocity. In contrast, emission from tidal material is mostly red- 
shifted. 

The depth of the H I absorption feature give an optical depth 
of 0.10 ±0.02 with a linewidth of 80 km s . Assuming a spin 
temperature of 100 K, this corresponds to an H I column density 
of (1.5 ±0.3) x 10 21 cm~ 2 . If the continuum source is located in 
the centre of the galaxy we would observe approximately half of 
the H I along the line of sight in absorption. This implies a total H I 

column density through the centre of the disk of (3.0 ±0.6) x 10 21 

—2 
cm . 



3.3 Molecular Gas 

Naturally weighted BIMA data were used to produce an image 
cube with a synthesised beam of Ofwhm = 8"8 x 5"3, correspond- 
ing to a physical scale of ~1.4 kpc. Figure [5] shows the integrated 
12 CO (J = 1—0) map overlayed on digitised sky survey S-band 
image. The CO data reveal four distinct peaks. The central peak is 
significantly brighter than the peaks observed offset to the east and 
west, its emission exceeding that in the eastern and western peaks 
by a ratio of 10 : 1 and 7 : 1 respectively. 

Whilst the molecular gas emission in NGC 5218 is clearly 
centrally condensed there is also a sizeable extended component. 
The central peak itself is extended in the easterly direction, well 
aligned with the optical bar axis. Examining the subsidiary peaks 
in the molecular gas, there is a clear peak located toward the east- 
ern end of the bar, whilst at the western end there is a larger area 
of emission extending from the bar end toward a region of star for- 
mation in the spiral arm. There is similarly a peak in the eastern 
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Figure 5. Integrated CO emission from BIMA 12 CO (J = 1-0) (0FWHM = 
8"81 x 5f'30 at a position angle of 13?0) observations of NGC 5218 over- 
layed on a digitised sky survey B-band image. The contour levels are (2, 4, 
8, 16, 32, 64, 128)Iybeam- 1 kms" 1 . 



region of the galaxy associated with a region of star formation in 
the eastern spiral arm. 

The integrated 12 CO (J = 1-0) flux for NGC 5218 was ob- 
tained from the primary beam corrected naturally weighted image 
cube. The molecular gas mass was estimated in solar masses from 
the integrated 12 CO (J = 1-0) flux using the empirical relation: 

M H , = 1.6 x 10 9 A 2 £> 2 J Sdv Mq 

where X is the wavelength in metres, D is the distance in mega- 
parsecs to the source, J Sdv is the integrated CO flux in Jy km s 
and CO-to-H 2 conversion factor, X = N(H 2 )/I C o = 2.8 x 10 20 
has been assumed (Bloe men et al ] |l986h . The 



observed integrated flux of 375 ± 75 J y km s i s in g ood agree 



K km s 

- J y 

ment with single-dish observations of Izhu et alj dl9991) . who ob- 
tained 363 ± 27 Jy km s , which implies that the BIMA observa- 
tions are not missing any flux. Our flux gives a total H2 mass of 
(6.9±1.4) x 10 9 M Q . 

The gas surface density in NGC 5218 is given by 
JShibatsuka et al.ll2003h : 



2h 2 



M pc~ 



:5.0x 10 2 cos(/) 



Ico 



Jy km s 1 arcsec 



and S gas = 1 .36Eh 2 . where i is the galaxy inclination a nd the facto r 
1.36 accounts for elements other than hydrogen (e.g., Alien ! 1973h . 
To estimate the peak central gas surface density from our 12 CO 
(J — 1—0) data we re-made the image cubes applying uniform 
weighting, giving a spatial resolution of 5"48 x 4"58. These data 
gave a peak central gas surface density of ~ 2000 Mq pc~ 2 . Ad- 
ditional o bservations o f the m olecular gas in this galaxy have been 
made by lOlsson etail d2004l) at higher resolution (3"75 x 3"54) 
using the OVRO millimetre array. Using these data they observe a 
peak gas surface density of ~ 3000 Mq pc~ 2 . 

The intensity weighted mean velocity field of the BIMA 12 CO 
(J — 1—0) data is shown in Fig. [6] These data are used later to fit 
a rotation curve for NGC 5218. The I2 CO (J = 2—1) emission ob- 
served with the JCMT is shown in Fig. [7] CO emission is detected 
over a region approximately 50" x 30". For comparison, we con- 
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Figure 6. Left: The intensity weighted mean velocity field of the BIMA 
12 CO (J = 1-0) data (Sfwhm = 8"81 x 5"30 at a position angle of 13.°0) 
for NGC 5218 The colour scale goes from 2716 (blue) to 3016 km s (red) 
and contour levels are (2730, 2750, 2770 . . . 3050) km s _1 . 
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Figure 7. Integrated 12 CO (J = 2-1) emission (Gfwhm = 20" 8) from 
JCMT observations of NGC 5218 overlayed on a digitised sky survey B- 
band image. The contour levels are (2, 6, 10, 14, 18, 22, 26, 30) K km s _1 . 



volved the 12 CO (J — 1—0) data down to the same resolution as the 
JCMT data (20" 8) and re-gridded the datasets to a common grid; 
the two datasets appeared consistent over their overlap region. 

We measure the 12 CO (J = 2- l)/ 12 CO (7=1-0) line ratio at 
the centre of NGC 5218 (RA 13 h 32 m 10?2, Dec +62° 46' 4"85) 
to be 0.8 ±0.2. This line ratio is similar to that found in Orion 
(Saka moto etaill 19941) and is consistent with the presence of gi- 
ant molecular clouds that have a kinetic temperature, 7^ ~ 20 K 
and densit y, njfy) ^ 10 3 cm 3 - T his ra ti° is similar to that ob- 
served by iBraine & Combesl ll992h in the centre of normal spi- 
ral galaxies (0.89 ± 0.06). The indication of a normal molecular 
gas environment at the centre of NGC 5218 is corrobo rated by the 
13 CO (J = 1 -0)/ 12 CO (J = 1 -0) line ratio of 9.0 ± 1 .3 dAalto et al.1 
which is comparable to the value observed in our Galaxy 
when averaged over the inner 300 pc. 



3.4 The relative gas content of NGC 5218 

The 12 CO (J = 1 -0) data for NGC 5218 give a molecular gas mass 
of (6.9 ± 1.4) x 10 9 M G (see SectionEOi. The mass of atomic gas 
still bound to the late-type system is (8.4 ±0.8) x 10 8 M , while 
the total atomic gas mass associated with the Arp 104 system is 
(2.8 ±0.2) x 10 9 M e (see Section^. 

Allowing for differences in the H2-to-CO conversion fac- 
tor, we compare the molecular-to-atomic gas ratio of NGC 5218 
to galaxies in the literature. Including all the atomic gas asso- 
ciated with the extended tidal features in NGC 5218, we ob- 
tain a molecular- to-atomic ratio, in terms of log(Me 2 /MHi) of 
0.40, which is significantly larger than that ob served in both nor- 
mal (-0.45 ± 0.07, from dCasoli et alj Il99fj) , and -0.15 from 



Bettonietal. 2003)) an peculiar SBb systems (-0.18 ±0.15, 



Casasola. Bettoni & G alletta 2004)). All mean ratios quoted in the 



literature indicate SBb systems have less molecular gas than atomic 
gas; in contrast, NGC 5218 appears to have at least twice as much 
molecular than atomic gas. 



3.5 The dynamics of the interaction 

The spatial extent of Arp 104's tidal tails can be used to make an 
approximate estimate the timescale of the interaction. The semi- 
circular tail extending to the north of NGC 5218 is ~ 6'.3, which 
corresponds to a physical distance of ~ 75 kpc, whereas the tidal 
bridge to the south is ~ 5', or approximately 60 kpc in extent. The 
two tidal extensions have comparable lengths, indicative of similar 
formation timescales. If we assume that the atomic gas associated 
with NGC 5218 was not accelerated during the encounter, then the 
velocities of the gas will be comparable to those prior to the in- 
teraction. The velocity of the atomic gas that has travelled approx- 
imately 70 kpc to the tip of the tidal tails can then be estimated 
from the rotational velocity of the system and the orbital velocity 
of the galaxy, v gas = (AV)/2±v, where v is the rotational veloc- 
ity of the galaxy from which the tidal tail originated and AV is the 
relative radial velocity of the two interacting galaxies. The interac- 
tion between NGC 5218 and NGC 5216 appears to be largely in the 
plane of the sky; the galaxies have very similar systemic velocities 
so we assume a relative radial velocity of zero. We estimate the ro- 
tational velocity of NGC 5218 to be ~ 247 km s (corrected for 
inclination). The spatial extent of the tidal tail and the estimated ve- 
locity of the gas imply a time since perigalacticon of approximately 
2.8 x 10 8 yr. 

As a further probe of the dynamical evolution of the Arp 104 
interaction we have modelled the observed H I distribution using 
collisionless Af-body simulations which provide a better constraint 
on the dynamics. Our observational data place some immediate 
constraints on the simulations, including: the absolute masses of 
the galaxies involved and their mass ratio; the relative velocity dif- 
ference of the line of sight velocities; the relative orientation of the 
galaxies; the physical separation of the galaxies and the length and 
orientation of the tidal tails. 

Dynamical, 6-band and A^-band mass estimates for NGC 5218 
and 5216, where appropriate, are given in Table[T] For the dynam- 
ical mass estimate of NGC 5218 a spherical mass distribution and 
inclination angle of 49° 5 (taken from the HyperLeda databas^H) 
were assumed. The 6-band mass estimates used a solar absolute 
6-band magnitude of 5.46 and assumed a mass-to-light ratio of 4 



see: http : //leda.univ-lyonl . It/ 
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Table 1. Galaxy Mass Estimates 



Galaxy 


fi-band 


if-band 


Dynamical 




Mass Estimate 


Mass Estimate 


Mass Estimate 




(M ) 


(Mq) 


(M ) 


NGC 5218 


8x 10 10 


7x 10 10 


1.3 x 10 11 


NGC 5216 


(4-7) x 10 10 


4x 10 10 





( Robe rts & Haynes| | 19941) for the late-type system and between 6 
ivan der Marellll99lh and 10 dFaber & Gallaghedll979h for early- 
type system. The A'-band mass estimates assume a stellar A'-band 
absolute magnitude of 3.32 and adopt a mass-to-light ra tio of 1 : 1 
for both the early and late-type systems dBell et alj2003h . Both the 
S-band and A'-band estimates are from the total emission magni- 
tudes of the galaxies. Whilst these mass estimates are subject to a 
large degree of uncertainty, they nonetheless provide a useful start- 
ing point for simulations of this system. 

Using these constraints a series of iV-body sim ulations were 
run using an in-house tree code ^Clemens et alJl99Sl) . Preliminary, 
restricted Af-body simulations were used to explore the parameter 
space. Simulations started with zero energy (parabolic) orbits. The 
restricted ,/V-body simulations allowed exploration of a range of 
possible perigalactic separations. Having established a parameter 
range within which simulations produced tidal features that were a 
good representation of those observed, further, self-consistent sim- 
ulations were run to confirm the validity of the selected parameters. 
The criteria for the 'best-fit' model was a simulated mass distribu- 
tion in which the projected separation of the galaxy pair was a good 
match to that observed and the morphology and velocity profile of 
the simulated galaxies most closely matched that of the H I distri- 
bution. 

The self-consistent A^-body simulations of galactic systems 
used two types of collisionless particles to simulate the stellar and 
dark matter; we ignore hydrodynamical effects in the formation of 
the tidal tails. Initial conditions for the simulate d late-type galax - 
ies are determined by a technique developed by Her nquisa J 1993h . 
Each late-type galaxy consists of a dark matter halo, a central stel- 
lar bulge, and a rotating disk. The particles are then positioned, 
following a spatial density profile modelled on observed profiles. 
The halo, disk and bulge components were assigned the mass ratio 
1.0: 1.0:0.3, the halo and bulge components we each composed of 
16384 particles, whilst the disk component contained 32768 parti- 
cles. The early -type system was simulated using the NEMO code 
( lTeubenlll995h to construct a Plummer model. In each case the 
model galaxies were evolved under gravity to confirm their stability 
prior to use in simulations. 

A series of simulations were run using models for an elliptical 
and spiral galaxy on zero energy (parabolic) or small positive en- 
ergy (hyperbolic) orbits. A number of mass scalings were explored. 
These included, adoption of the A*-band mass estimates, adoption 
of the A'-band mass ratio but using the mass scale defined by the dy- 
namical mass estimate for NGC 5218 and adoption of the dynami- 
cal mass estimate for NGC 5218 and the A'-band mass estimate for 
NGC 5216. The last of these combinations implied a mass ratio of 
the late to early-type galaxies of 4 : 1 . These simulations produced 
short tidal tails, unlike those observed in the H I data. On the basis 
of this work we suggest that the mass ratio of the two galaxies is 
closer to 2 : 1, as implied by the A'-band mass ratio. 

All the simulations run using a 2 : 1 mass ratio for the late to 
early-type galaxies produced a tidal-tail morphology similar to that 



observed in the H I data. The results obtained appeared relatively 
robust. Several simulations were run, with an initial disk inclina- 
tion of 0, 30, 60 and 90 degrees. The best-fit set of parameters 
were based on a parabolic interaction with a perigalactic separa- 
tion of 12.5 kpc, using galaxies with a masses approximated using 
the A'-band mass-to-light ratio; the simulated data, for an initial 
disk inclination of 60 degrees, are shown in Fig. [8] The time from 
perigalacticon to the point at which the simulated data best repre- 
sented the H I observations was ~ 50 time-steps, corresponding to 
a timescale of ~ 3.3 x 10 s yr. This is in good agreement with the 
crude dynamical estimate derived above. 



4 ANALYSIS OF THE BAR IN NGC 5218 
4.1 Rotation Curve 

The first step in the analysis of the bar is to attain a rotation curve 
for NGC 5218. The rotation curve was established predominately 
using the VLA H I data. Using the channel maps (Fig. O, spatial 
and velocity selection were use to isolate gas which is dynamically 
associated with NGC 5218, this corresponded to gas out to a radius 
of approximately 9 kpc. NGC 5218 contains a bright central radio 
continuum source and as a result H I in the nuclear region is seen 
in absorption, making it difficult to fit a rotation curve to the H I 
data in this region. However, in the central region our 12 CO (J = 
1—0) observations provide a reliable estimate of the rotation curve. 
We therefore use H I data for radii beyond about 2.5 kpc and CO 
data within this radius. The 12 CO (J = 1—0) observations have the 
additional advantage that they are higher resolution than the H I 
data, re ducing the proble ms of beam-smearing in the centre of the 
galaxy iSofue et al.ll 19991) . 

Rotation curves for NGC 5218 were derived from the 
intensity-weighted velocity field of the H I and CO data using 
the AIPS task GAL. For both the H I and CO datasets we used 
the highest resolution images, so as to minimise the effect of 
beam smearing. For the Hi data this corresponded to a spatial 
resolution 13"3 x 10"7, and for the CO data the resolution was 
8"8 x 5"3. The rotation c urves were derived using a tilted-ring 
model JWarneretal.ll 19731) . This method approximates the galaxy 
as a set of concentric rings of gas at increasing radii, characterised 
by a circular velocity, position angle and inclination angle. At- 
tempts to derive the rotation curve leaving the dynamic centre and 
inclination angle as free parameters yielded unphysical results. The 
dynamical centre was therefore calculated by performing a global 
rotation curve fit to Hi datasets at three different resolutions for 
both Brandt and Exponential rotation curves independently. The 
dynamical centre was found to be stable across all fits and was 
therefore fixed for the purposes of the tilted-ring curve fitting. Op- 
tical data were used to fix the inclination angle which was taken 
from the Hyperleda databasfl 

The H I rotation curve was fitted by breaking the galaxy into 
annuli 5" in width (approximately the half-beam width) and fitting 
the velocity field in each annulus, keeping the dynamical centre and 
inclination angle fixed. Radii below 10" in which the H I is seen 
predominately in absorption were excluded. A similar fit was per- 
formed on the molecular gas data, using the same fixed dynamical 
centre and inclination angle. For these higher resolution data 2"5 
annuli were used, and the data were fit from zero radius outwards. 
In Fig.|9]we plot the circular velocity v(R), resulting from analysis. 

3 See http: //leda.univ-lyonl .f r/ 
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Figure 8. The distribution of collisionless particles obtained from a simulation of the encounter between NGC 5216 and NGC 5218. Red particles indicate the 
disk particles of the late-type component and green particles indicate the particles comprising the early-type component. The orbital plane of the galaxies is in 
the plane of the figure. The axes indicate distance in kpc. 
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Figure 9. The model circular velocity de-projected on the galaxy plane 
as obtained by fitting a tilted-ring model to the combined molecular (plus 
signs) and atomic (diagonal crosses) gas data for NGC 5218. 



Figure 10. The angular velocity curve Cl(R) (stars) and derived curves Q. + 
k/2 (diagonal crosses) and £1— k/2 (plus signs), using the model circular 
velocity from the combined molecular and atomic gas datasets. 



We fit a function to these data the differential of which was evalu- 
ated to allow computation of K, the epicyclic frequency. In Fig. 1101 
we plot the angular velocity, Q.(R) = v(R) /R and £1 ± k/2 curves. 
The error bars on Figs [9] and [10] do not include any possible bias 
due to beam smearing. 
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4.2 The length of the bar and co-rotation 

Science-grade images in J, H and K s were downloaded from the 
Two Micron All Sky Survey (2MASS) archivtQ. The J, H and 
A" s passbands are less sensitive to extinction by dust in the cen- 
tral regions of the galaxies than shorter wavelengths, and offer a 
bar tracer that is unlikely to be contaminated by recent s tar forma- 



tion. We adopt the method oflRegan & Elmegreenldl997 ) - see also 
iMartini et al] J200lh : [sheth et aljj2002t) and lSheth et alU2003l) for 



examples of this technique - to estimate the length of the bar; this 
assumes the bar ends correspond to the point of maximum elliptic- 
ity when fitting the galaxy isophotes. The ellipticity of the 2MASS 
image isophotes, in the plane of the sky, were fit using the IRAF 
STSDAS task ELLIPSE, which measures the surface brightness, el- 
lipticity and position angle of the semi-major axis. FigurefTTIshows 
the ELLIPSE fitted ellipticities for both the J- and AT-band 2MASS 
images. Both images give very similar results; the /-band image in- 
dicates a maximum ellipticity of 0.618 ± 0.008 at a radius of 13" 3 
and the AT-band image a maximum ellipticity of 0.610 ±0.008 at a 
radius of 14"0. 

The bar's de-projected length / can then be obtained from 
the observed length / n h s using / t, s cos 2 4> + sin <5cos 2 i 

dAmaboldiet"aIlll995h . where i is the inclination angle of the 
galaxy (49.°5) and <I> is the angle the bar makes with galaxy's ma- 
jor axis in the plane of the galaxy. To determine S> we take the 
position angle of the galaxy to be the average of position angles 
fitted to the H I data using the tilted-ring model in GAL and the 
position angle of the bar to be that obtained by the IRAF package 
ELLIPSE at maximum ellipticity. These angles are 71°4±4°3 and 
81°6±0°3 respectively. Thus, the angle the bar makes with the 
major axis in the plane of the sky <5' = 10° ± 5°, and <I> follows 
since tan(<&') = tan(<J>)cos(i), thus $ = 15° ±5°. These angles 
yield a very small correction factor in converting between the ob- 
served and de-projected bar length: / = 1.02 x l ^ s , gi ving / = 28"6 
and r = 14 "3. Assuming co-rotation is at 1.2 x R iAthanassoulal 
Il992lll994l) yields a co-rotation radius of 17" 1 which corresponds 
to a physical radius of 3.4 kpc and gives and estimate of the pattern 
speed of the bar of 40 k m s -1 kpc -1 , which is simila r to that seen in 
other SB galaxies (e.g. , iMe rrifield & Kuijke^ 19951; jGerssen et all 



1 19991 ; lLaine et al.|[l998h and in our Galaxy (i Dehnenll 19991) . Exam- 
ining the angular velocity obtained from the rotation curve fitting at 
a radius of 3.4 kpc we find the velocity is significantly larger than 
the peak value of the Q. — k/2 curve; this implies there is no inner 
Lindblad resonance in NGC 5218. We note that the calculated ra- 
dius of co-rotation (3.4 kpc) is not dissimilar to the radius obtained 
by measuring distance between the molecular gas peaks toward the 
ends of the bar (3.9 kpc). 

Constraints on the existence or not of an inner Lindblad res- 
onance in NGC 5218 depend sensitively on Q.(R) and k(R). As 
such they are subject to a large degree of uncertainty. In particular, 
the absence of an observed peak in the £2 — k/2 is symptomatic 
of a mass distribution which is not strongly centrally condensed; 
inner Lindblad resonances only develop if the mass distribution is 
sufficiently centrally peaked. It is well known that low resolution 
observations smooth out the steep central rise of rot ation curves 
which results from the core and bulge of the galaxy 1 Sof ue et"al1 
1999). It is possible that the limited resolution of our observations 
is making it difficult to detect the steeply rising rotation curve in 



4 Available at 

http : //irsa. ipac . caltech.edu/applications/2MASS/LGA/ 
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Figure 11. The ellipticity of NGC 5218 as a function of radius obtained 
from the ellipse fitting of the IRAF STSDAS task ELLIPSE, for both the 
2MASS /-band (green) and K-bwd (red) images. 



the central region of NGC 5218 and as such yielding an artificially 
smooth £2 — k/2 curve. 



4.3 Comparison of the H I and CO dynamics 

Figure|9]shows the fitted circular velocities for both the H I and CO 
datasets. Around radii of approximately 3 kpc, where we have both 
atomic and molecular gas observations there is a clear deviation 
in the velocity structure of the two components. This is not unex- 
pected. The fall in the angular velocity of the molecular gas around 
3 kpc, in the region of the bar end may be a result of streaming mo- 
tion s within the bar. Gas on the elongated, so-called x\ orbits in the 
bar l lAthanassoula 1994), has a rotational velocity at the apex of the 
ellipse, near the end of the bar, which is lower than that of gas on a 
circular orbits at the same radius, therefore its £l(R) will be lower. 
By comparison, if the H I has a large scale height than the molec- 
ular gas, it is likely to have been less influenced by the bar in this 
region. We have also performed a more detailed comparison of the 
Hi and CO dynamics by comparing the Hi and 12 CO (J = 1—0) 
spectra at four locations across the disk. The spectra were extracted 
from data cubes at a common resolution of 17"2 x 13"5. In the 
central region the linewidths and lineshape of the spectra were very 
similar. In the more extended regions the molecular gas displays a 
broader linewidth than the H I. In particular, at the eastern end of 
the bar, the molecular gas displays a two-peaked profile with the 
atomic gas indicating emission in only one of these peaks. 



5 DISCUSSION AND CONCLUSIONS 
5.1 Enhanced central molecular gas 

In Section [3~l4l we showed that NGC 5218 has a large molecu- 
lar gas mass for an SBb galaxy. The molecular gas distribution is 
strongly centrally concentrated with a peak surface density greater 
than 2000 Mp, pc~ 2 (Section l3.3t . Numerous numerical simula- 
tions predict that bars, whether spontaneously or externally trig- 
gered, can induce large scale gas inflow toward the central region s 
of galaxies (e.g., MH96; lNoguchill9 88; Barn es & Hernquistll99ll) . 
Two recent obser vational studies have fo und strong evidence for 
this process; both lSakamoto et alj dl999al) and lSheth etal] fc005l) 
have demonstrated that barred galaxies have a larger concentration 
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of molecular gas in thei r central kiloparsec than unbarred systems. 
ISakamoto etaU dl999tj) used the NRO-OV RO survey of 20 nearby 
spiral galaxies (see lSakamoto et al . 1999b) to show that the degree 
of concentration of molecular gas in the central kiloparsec of barred 
spiral galaxies is greater than in their unbarred counterparts. They 
quantify this effect using the ratio of the molecular gas surface den- 
sity in the central kiloparsec to that over the entire optical disk 
(/ con ), finding this ratio is on average four t imes larger in barred 
than in unbarred systems. ISheth et all d2005h studied a sample of 
50 nearby spirals - including both late type barred and unbarred 
spirals, and CO-bright and CO-faint galaxies - and found the effect 
of the bar equally pronounced in late and early Hubble type spi- 
rals. Atomic gas observ ations also indicate higher gas colum n den- 
sities in barred galaxies. iRisaliti. Maiolino & Salvatil lfl999) exam- 
ined the atomic gas column density in front of Seyfert nuclei using 
X-ray measurements and found significantly higher column densi- 
ties in Seyferts with barred hosts than those in the unbarred sam- 
ple. To estimate / con for NGC 5218 we remade the channel maps at 
the highest possible spatial resolution (6fwhm ~ 5", ~ 1 kpc). Al- 
though beam smearing introduces uncertainty into our analysis, we 
find a 12 CO (J = 1 -0) flux over the central kiloparsec in NGC 5218 
which corresponds to a gas mass of ~ 1.5 x 10 9 M . This is equiv- 
alent to a molecular hydrogen surface density of ~ 1300 Mq pc~ 2 
and a gas surface density of ~ 1800 M Q pc~ 2 . Thi s is approx- 
imatel y six times the mean of the barred sample of She th et all 
J2005h (~ 300 M© pc~ 2 ). In addition, the data for NGC 5218 in- 
dicate an / con value of 113, almost twice the average of the BIMA 
sample (~ 60) and slightly larger than the average of the NRO- 
OVRO sample (~ 100). The discrepancy between the two litera- 
ture values can be attributed to the NRO-OVRO sample consisting 
mainly of CO-luminous galaxies. In contrast, the BIMA sample 
was selected on the basis of B-band magnitude and thus represents, 
at least from the perspective of CO flux comparison, a less biased 
sample. The existence of the bar in NGC 5218 has led to a very 
large central molecular gas surface density, even within the context 
of barred galaxies. For comparison we note that typical mean den- 
sities of the largest molecular cloud complexes in M31, M33 and 
M5 1 are in the range 40—500 Mq pc~ 2 , whilst the gas surface den- 
sity of a typical giant molecular cloud in the oute r disk of the Milky 
Way is 200 M pc -2 dScoville & Sandersll 19871) . 

The impact of the bar potential on the gas distribution in 
NGC 5218 is consistent with the bar's ellipticity, which implies it 
is a strong bar. The ellipticit y of a bar is th ought to provide a good 
measure of the bar strength. iMartinl (T995) identify strong bars as 
those with an ellipticity greater than 0.4; the 2MASS A"-band data 
for NGC 5218 give a deprojecte d ellipticity for NGC 52 18 of the 
bar of at least 0.4. Simulations bv lRegan & Teubenl d2004l) find that 
weak bars have almost no effect on the radial distribution of the 
ISM. In contrast, strong extended bars are likely to be more effec- 
tive at driving gas toward the centre of the galaxy. Their simulations 
also indicate that bars are only able to drive gas inwards to a radius 
at which a ring forms; for bars containing inner rings, little net in- 
flow takes place inside this radius. Nuclear rings are likely to form 
when so-called xi orbits exist; such orbits require the presence of 
two inner Lindblad resonances. Our very high central gas densities 
are consistent with no inner ring and no inner Lindblad resonances 
as our analysis of the bar dynamics suggests. 

5.2 CO-to-H2 conversion factor 

The molecular gas mass in NGC 5218 and the resulting gas col- 
umn densities are calculated using the CO-to-H2 conversion fac- 



tor^ = AT(H 2 )//co = 2.8 x 10 20 cm" 2 K km s" 1 dBloemen et all 
1986), which is appropri ate for molecular gas clouds in the Milky 
Way dSanders et alj 1984b . In adopting this conversion factor we as- 
sume that the conditions of the molecular gas clouds in NGC 5218 
are similar to the molecular gas clouds in the Milky Way. This is 
an assumption, and it is possible that the CO emissivity in this sys- 
tem is enhanced due to high pressur es and/or velocity dispersion in 
the bulge region dReean et alj200ll) . However, the obse rvation of a 
13 CO (J= 1-0)/ 12 CO(7= 1-0) line ratio of 9.0± 1.3 dXalto et al.1 
Il99ll) in NGC 5218, which is comparable to the value observed in 
the Milky Way when averaged over 300 pc, lends weight to the 
adoption of the Milky Way conversion factor. 

Although the molecular gas distribution in NGC 5218 is cen- 
trally concentrated, there is also evidence of more extended emis- 
sion (see Fig. [5)- In particular, the 12 CO (J = 1—0) integrated emis- 
sion indicates two molecular gas condensations slightly offset from 
the bar ends. The radius of these gas peaks (3.9 kpc) is similar to 
the calculated radius of co-rotation (3.4 kpc), see Section |4~2l The 
presence of these gas complexes, like that of the large central peak, 
can be attributed to the existence of a barred potential. Gas con- 
centr ations at bar ends have been observed in a number of galaxies 
(e.g. iKennev & Lordfl99ll; bownes et all 19961 : ISheth et alj|200d 
120021) . These concentrations may be stagnation points in the gas 
flow; gas on x\ orbits at the bar ends will have a smaller angular 
velocity than more central gas due to its larger radius. Their pres- 
ence may also be due to the convergence of gas on elliptical orbits 
in the b ar with gas on more cir cular orbits in the region just beyond 
the bar dKennev & Lord 1991); these gas complexes are located in 
a transition region between the spiral arms and bar ends. Although 
the bar in NGC 5218 appears to have contributed to a central con- 
centration of gas, this process appears to be ongoing. 

We now consider the origin of NGC 5218's high molecular 
gas mass. In Section |3^4l we showed that NGC 5218 has a compar- 
atively low H I mass. The bar may have led predominately to the 
transport of existing molecular material from the disk to the centre 
of NGC 5218, or this transport could have involved both atomic 
and molecular gas with subseq uent conversion of atomic t o molec- 
ular gas in the central region. iDavoust & Contiriil d2004l) studied 
the effects of bars on the neutral hydrogen content of starburst and 
Seyfert galaxies and found that the barred galaxies had a lower to- 
tal H I mass than unbarred galaxies. They suggest that atomic gas 
is funnelled toward the centre of barred galaxies where it is then 
converted to molecular gas. In Section [3T2I we used H I absorption 
observations to estimate the atomic gas column density in the centre 
of NGC 5218 to be (3.0 ±0.3) x 10 21 cnr 2 . Various studies indi- 
cate a 'transition' column density at which H I is sufficiently dense 
for self-shielding to allow the formation of H2 of ~ 5 x 10 20 cm~ 2 , 
increasing to ^ 10 21 cm~ 2 for large molecul a r clouds dSavage et al.l 



19771; iFederman. Glassgold & KwarJ I l979l ; iFranco & Coj|l986t 



Reach. Koo & Heilesl ll994). The H I column density observed at 
the centre of NGC 5218 lies above these values so H I to H2 con- 
version may be taking place at the centre of NGC 5218. The high 
H2 to H I ratio in the centre of NGC 5218 and the dynamical consis- 
tency of the two components in the central region (see Section |4~3l l 
is certainly consistent with inflow of both components and efficient 
H I to H2 conversion. 

5.3 Star formation 

In Section |3T| we estimated the SFR of NGC 5218 using both radio 
and FIR data; these data yield SFRs of 3.4 and 5.3 M p yr~' re- 
spectively. A large proportion of the star formation is taking place 
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in the dusty central region of NGC 5218 (see below) where we ob- 
serve a very large molecular gas density. However, both the radio 
and FIR provide largely extinction-free probes of the star formation 
and as such should yield robust estimates of the SFR. Comparison 
of the NVSS (9fwhm = 45") flux densit y (30.4 ± 1 .0 mjy) and the 
FIRST teecker. White & Helfan dll 19941) (0fwhm = 5 "4) flux den- 
sity (23.06 ±0.14 mjy) suggests at least 75% of the galaxy's star 
formation is coming from the central kiloparsec region. 

Before discussing the star-formation properties of NGC 5218 
we note a potential caveat with the estimation of the SFR of this 
system. A number of authors have examined the optical spec- 
troscopy o f NGC 5218 for AGN activi t y. O n the basis of their in- 
vestigation Ash bv. Houck & Matthews! jl995h conclude the emis- 



sion li ne properties of NGC 5218 are ambiguous, and lVeilleux et al.l 
dl995h classify it as a LINER. NGC 5218 sits com fortably on the 
radio-FIR correlation dYun. Reddv & Condonll200ll) with no obvi- 
ous radio or FIR excess; the two SFR estimators yiel d very sim- 
ilar re sults (Section [3. It . In addition, recent work by iDudik et al.l 
(2005) suggests that in their sample of IR-bright LINERS the low 
X-ray luminosities and Eddington ratios imply the contribution 
from the AGN to the FIR luminosity is likely to be small. How- 
ever, we cannot exclude the possibility that nuclear activity in this 
system leads to an overestimated FIR SFR. 

5.3.1 Is the level of star formation large? 

NGC 5218 appears to have a substantial star formation rate. 
The 60-^lm luminosity of NGC 5218, 1.5 x 10 1() L , is more 
than three times the med ian for Sb-Sbc galaxies found by 
iDevereux & Hameedl d 19971) . although it does not qualify it as Lu- 
minous Infrared Galaxy (LIRG). NGC 5218's //-band normalised 
SFR, an approximation o f star formation per unit m ass, implies 
a slightly enhanced value dDever eux & Hameed 1997 ). Th e mean 
value of log(LpiR/LB) obtained by Devereux & Hameedl dl997h 
for Sb-Sbc galaxies was —0.48 ± 0.34, whereas the value for 
NGC 5218 is —0.05, consistent with a high current SFR in relation 
to its recent star formation history (although it may also reflect, to 
some extent, increased gas and therefore dust in its centre). 

The star-formation efficiency (SFE) is the ratio of the current 
star-formation rate to the mass of molecular gas: for NGC 5218 the 
SFE derived from the FIR star formation rate and total molecular 
gas mass is ~ 8 x 10~ 1Q yr -1 . This value, whilst larger than the 
average for Sb galaxies (Murgia et al. 2002), is still si gnificantly 
below the ~ 30% over 10 8 yr found in starburst systems dKennicutil 
Despite NGC 5218's relatively high star-formation rate, its 



global star-formation efficiency appears to be more in keeping with 
that of a normal Sb galaxy than a starburst system. The star forma- 
tion efficiency of the central 2—3 kpc is approximately three times 
the global value, falling to ward the lower end of those systems cla s- 
sifie d as nuclear starb ursts Planes as. Colina & Perez- 01eaUl997l) 

ISakamoto et ail dl999al) have investigated the gas-to- 
dynamical mass ratio in the central kiloparsec of a sample of 20 
nearby spiral galaxies. They observe a larger ratio in galaxies with 
an H II spectral classification, and interpret this as evidence for 
more active star formation and nuclear spectra dominated by H II 
regions in galaxies which have large amounts of gas given their 
size. They suggest that the effect may be related to the stability of 
jas di sk, relating Mg as /Md yn to the Toomre 



the nuclear molecular 



Q parameter dToomrd 19641) . For comparison with these data we 
use the highest resolution CO data set for NGC 5218 to examine 
the gas-to-dynamical mass ratio in central kiloparsec. We integrate 
the flux obtained in the central molecular gas peak and compare 



the mass obtained with an estimate of the dynamical mass in this 
region. We select a physical region with a radiu s of r ~ 1 kpc, 
larger t han the region (r = 500 pc) examined by ISakamoto et al.l 
d 1999a ) because we can make a more robust estimate of Mgas/M^y,, 
on this scale given our available data. The region examined lies 
well within the physical extent of the bar; the molecular gas 
dynamics in this region are dominated by circular motion (Fig.[6}. 
Assuming Keplerian rotation and a spherical mass distribution, the 
dynamical mass is 

\ M J \kpc ) Vkms- 1 

where v is the rotational velocity in the plane of the disk at a radius 
r. The 20% linewidth over the 1.1 kpc radius region is 280 km s . 
Correcting for the inclination of the galaxy yields a rotational ve- 
locity in the plane of the disk of 184 km s -1 , corresponding to a 
dynamical mass of 8.5 x 10 9 Mq. Integrating the CO flux observed 
over this area yields a molecular hydrogen mass of 1.6 x 10 9 M 
and a molecular gas mass of 2.2 x 10 9 M s . The molecular gas mass 
fraction over the central 1 . 1 kpc r adius region in NGC 521 8 is 26%, 
lo" above the value obtained by I Sakamoto et al] dl999al) for H II 
dominated systems (17.6 ± 8.0% over the central kiloparsec). If we 
reduce the physical scale over which we undertake the comparison, 
the mass fraction increases dramatically. Our data imply that the 
gas-to-dynamical mass ratio over the central region of NGC 5218 
is large when compared with actively star forming H II systems in 
the Sakamoto et al. sample. 

We note that the derived Mg as /Miy n is only related to the 
Toomre Q parameter and the central gas stability if the gas is in 
the form of a nuclear disk. Bar-streaming motions introduce in- 
herent uncertainty into the use of the ratio Mo^/M^ yn as a probe 
of the stability. In addition, many theorists are sceptical about us- 
ing the Toomre Q parameter as a measure of instability toward star 
formation, arguing that feedback from such st ar formation w ould 
low er the parameter o n short timescales (e.g., |Combes 20 0j], see 
also Koda et al. 2005). Nonetheless, the results of Sakamoto et al. 
( 1999a) indicate that all except one of the H II galaxies in their sam- 
ple have an M gas /M(jy n greater than 0.1 and we believe comparison 
with these observations is useful. 



5.4 Bar formation in NGC 5218 

Both simulation-based and observational evidence exists linking 
tidal perturbations with the triggering of bars. In particular, ob- 
servations in the literature point to bars in early-type galaxies, of 
which NGC 5218 is an example, being produced in this way. This 
result is supported by simulations which indicate that a bulge com- 
ponent acts to stabilise the galaxy against spontaneous bar for- 
mation, an increasingly massive bulge offering increased stability 
(MH96). It may be that for galaxies with bulges an external trig- 
ger is a prerequisite for a bar to form. To investigate any link be- 
tween the bar in NGC 5218 and its tidal interaction, we look in 
more detail at the nature of the bar, undertaking comparisons with 
both simulation-based and observational data. 

Simulations suggest that interactions can affect bar formation 
in galaxies in a number of different ways: they can trigger bar for- 
mation in an otherwise stable galaxy; speed up the process of ba r 
formation in a system which is already u nstable dGerin et alj 1990h ; 
regenerate a dissolved or weakened bar dBerentzen et a l. 2004) and 
enhance or decrease the strength of an existing bar ICjerin et all 
1990). All such scenarios are possible in NGC 5218, although the 
strength of the bar, as measured from the axial ratio, would imply 
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that the tidal encounter is unlikely to have significantly weakened 
the bar. 

In an investigation of th e blue and near-infrared surface pho- 
tometry of 15 barred spirals. lElmegreen & Elmegreenl dl985h ob- 
served two different types of major axis surface density profiles, 
expon ential and flat. On the basis of their simulations, iNoguchil 
argue that exponential profiles are found in spontaneously 
produced bars, whilst flatter profiles occur in bars produced by 
tidal instability. Unfortunately the photometric data available for 
NGC 5218 are not of sufficiently high resolution to confidently 
discern between the two profiles; the large central peak observed 
in all bars yields a profile that appears exponential when smoothed 
to lower resolution. Nonetheless, the SDSS z-band light profile in- 
dicates some flattening and a flattened profile is certainly possible. 
In Noguchi's interpretation, this implies the bar profile is consis- 
tent with a bar produ ced by tidal perturbation. Recent work by 
iBerentzen et al.l J2004) has questioned this in t erpret ation, arguing 
that the results of lAthanassoul a~& Misiriotisl fe002h indicate that 
both profiles occur in isolated systems and that the different profiles 
are associated with varying bar strengths rather than the triggering 
mechanism, flatter profiles being produced in stronger bars. 

Our analysis of NGC 5218 gave an e stimated pattern speed for 
the bar of ~ 40 km s~ 1 kpc~ 1 . Studies bv lMiwa & Noguchil Jl998h 
show that in weak tidal perturbations the properties of the bar de- 
pend largely on the internal structure of the host ga laxy, rather 
than t he speed of the perturber (see also simulations bv lGerin et al.l 
1990). However, for stronger tidal perturbat ions the pattern spee d 
is dependent on th e perturber mass (see also Berentzen et alj2004h . 
iMiwa & Noguchil argue that for the strongest perturbations tidally 
induced bars can rotate slowly with an inner Lindblad resonance 
near the bar ends, whereas bars produced by weaker perturba- 
tions have no inner-Lindblad resonances, rotate more rapidly, and 
end near co-rotation. In the case of NGC 5218 we inferred that 
there are no inner-Lindblad resonances, and that if the bar-ends are 
co-rotating its angular velocity is ~ 40 km s kpc -1 . However, 
since the £2 — k/2 curve provides a relatively constant resonance 
throughout the galactic radii at a velocity of 10 km s kpc~' , we 
cannot exclude the possibility that the bar in NGC 5218 is rotating 
at this, much slower velocity, and ending before co-rotation. 

A number of investigations of tidal triggering of bar forma- 
tion have looked at the relationship between strength of interactio n 
and bar formation [Mi wa&Noguchi| [l998; Beren tzen etai1l2004l) . 
Quantitative comparison between these simulations and our results 
is difficult since many of the results are sensitive to the mass con- 
centration of the disk and the disk-to-halo mass fraction which are 
difficult to constrain observationally. Nonetheless, we can make 
some important observations about this system. Given the appar- 
ently low mass ratio of NGC 5216 to NGC 5218 (1 : 2 using A"-band 
data, Section [3~5l >. the perigalactic separation required to produce 
a tidal perturbation of sufficient strength to induce bar formation is 
likely to be small. The necessity of a small perigalactic separation 
is consistent with the simulations, which indicate that production of 
the extended tidal features observed in the H I distribution require 
a small perigalactic separation. 

The gas surface density in the central kiloparsec of NGC 5218 
(1300 Mp, pc~ 2 ) is significantly larger than the mean of both the 
unbarred (107 ±29 M B pc~ 2 ) and barred sample (309 ±71 M Q 
pc~ 2 ) of Sheth (2001). The enhanced surface density observed in 
NGC 5218 can be accounted for by an additional molecular gas 
mass of the order of 9 x 10 8 Mq. Estimates of gas inflow rates 
in ba r red galaxies vary wi d ely (0.1 — 17 M w yr^.lFriedri & Benzl 
1 19931 : lOuillen et alj 1 19951 ; lLaine et al.l 1 19981 ; iRegan et al.l 1 19971 ; 



ISakamoto et al.l 1 1999al; llono et alj liooj) ; adopting the value in- 
ferred by lOuillen et alj h995h of 4 ± 2 M© yr 1 via calculation 
of the torques on the molecular gas using near-IR images, we find 
the time taken to produce the observed enhanced central gas density 
in NGC 5218 is similar to the time since perigalacticon. However, 
it must be remembered the inflow rate is not certain, and may vary 
with time. 



5.5 Conclusions 

NGC 5218 is a barred interacting spiral galaxy which has under- 
gone a prograde encounter with an elliptical, NGC 5216, sometime 
in the recent past. Both the H I velocity field and the simulations in- 
dicate the timescale since perigalacticon is of the order of 3 x 10 8 
yr. NGC 5218 has a very large central molecular gas surface density 
when compared to the molecular gas surface density i n the centra l 
kiloparsec of the BIMA sample of unbarred galaxies dShethll200ll) . 
The difference can be accounted for if we assume gas inflow to the 
central region over the time since perigalactic on and at a rate com - 
parable to those inferred in other systems (e.g., lOuillen etal J 19951) . 
Our observations are consistent with the interaction of NGC 5218 
with NGC 5216 triggering an instability in the galactic potential or 
enhancing an existing instability that has led to gas inflow yielding 
a central molecular gas density at the elevated level observed. 

The nature of the bar, molecular gas and dynamical time scales 
we infer for NGC 5218 are very similar t o those in UGC 2855, as 
discussed by Htit temeister, Aalto & Walll dl999h . This interacting 
galaxy also has a high central gas concentration with no evidence 
for any inner Lindblad resonance or shocking of the molecular gas. 
Htittemeister et al. deduce that the bar is very young and the system 
is in a pre-starburst phase. 

NGC 52 1 8 has a relatively large SFR (Section|53ll It has a nor- 
mal disk-wide star-formation efficiency although the nuclear region 
has a star-formation efficiency toward the lower end of starburst nu- 
clei. It is possible that NGC 5218 is in the early stages of bar-driven 
evolution and has not yet reached its period of most active star for- 
mation. Whilst the bar in NGC 5218 appears to have contributed to 
large central molecular gas density there is still a significant frac- 
tion of the molecular material found in the more extended regions 
of this system. The process of gas transfer toward the central re- 
gion of NGC 5218 is clearly ongoing. The molecular gas density in 
the centre of NGC 5218 may not yet be sufficiently large to have 
triggered a starburst phase. We suggest that NGC 5218 is likely 
to undergo a period of more enhanced star formation in the near 
future. This proposition is supported by the following evidence. 

• The large concentration of molecular gas and very high sur- 
face density in the central region of NGC 5218. 

• A sizable, extended molecular gas component which indicates 
that the process of gas transfer to the central region is ongoing. 

• NGC 5218's large FIR-luminosity and SFR when compared 
with other galaxies of the same Hubble type. 

• The large Lm/Lg ratio, indicating that the current-to-recent 
star formation in NGC 5218 is larger than an average early-type 
spiral. 

• The large M sas /M^y n ratio in the centre of NGC 5218 imply- 
ing a large quantity of gas in this region given its size. This may be 
indicative of instability toward star formation. 
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